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Abstract 

The concept of point interactions offers one of the most suitable guides towards 
a quantitative analysis of properties of certain specific non-Hermitian (so called 
VT— symmetric) quantum-mechanical systems. This is illustrated on a double-well 
model, an easy solvability of which is shown to lead to a clear picture of the mecha- 
nisms of the unavoided level crossing and of the spontaneous VT symmetry breaking 
at a certain strong-non-Hermiticity boundary. Below this limit the model is shown 
suitable for an explicit illustration of technicalities related to the Standard probabilis- 
tic physical interpretation of bound states in VT — symmetric quantum mechanics. 
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1 Introduction 



A consistent probabilistic interpretation of bound states in quantum mechanics 
requires that their norm is conserved in time [1]. In the light of the well known 
Stone's theorem this means that their time-evolution must be unitary so that the 
generator of this evolution (= the physical Hamiltonian H p ) cannot be, apparently, 
non-Hermitian. Still, the concept of a " non-Hermitian Hamiltonian" H n ^ is in a 
more or less current use in many mathematical considerat ions (e.g., in perturbative 
calculations [2] or in the schematic models of gravity [3] or supersymmetry [4]) as 
well as in the various phenomenological studies (e.g., in nuclear physics [5], solid- 
state physics [6], particle physics and field theory [7] or in the model-building in 
magnetohydrodynamics [8] and cosmology [9]). 

In such a context the use of the solvable point interactions seems to offer a 
natural background for a simplification of some mathematical challenges [10] or of 
technical problems [11] as well as for an improvement of the feasibility of numerical 
calculations [12] and for a deepening of the related physical interpretation of the 
underlying dynamical processes [13, 14]. In our paper we shall address several aspects 
of non-Hermiticity in this framework, especially those which were recently studied 
using the comparatively less easily tractable piecewise constant potentials [15]. We 
feel encouraged by the success and impact of these studies which already found 
applications in supersymmetric quantum mechanics [16] and which seem to exhibit 
certain promising phenomenological (e.g., localization) properties [17]. 

We are persuaded that the class of the point-interaction models could be even 
better capable to throw new light on some of the paradoxes connected with the 
non-Hermiticity. We are going to pay attention to one of the most elementary rep- 
resentatives of this class which is just a simplified form of the potential which we 
already analyzed in our older, purely numerical study [12] (cf. also [13]). As long as 
all the similar models combine solvability (i.e., simplicity) with flexibility (controlled 
by several real parameters) we believe that there are good chances of an explicit 
confirmation of their full compatibility with the Standard quantum theory. Here we 
are going to perform the first steps in this direction. 
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2 A Zoo of non-Hermitian Hamiltonians 



Our present study is motivated by the fact that in the light of our introductory 
comment the use of non-Hermitian Hamiltonians sounds like a paradox. Still, in the 
literature one may notice an emergence of various sophisticated possibilities of mak- 
ing many of these models compatible with the concepts of the physical observability 
[5, 18], with the Standard principle of correspondence [19, 20] or with an introduction 
of relativistic kinematics [21] etc. 



2.1 Quasi-Hermitian Hamiltonians 

In a preliminary remark let us remind the reader that even many finite-dimensional 
non-Hermitian matrices A ^ may be re-interpreted as Hermitian after one re- 

— * 

defines the scalar product of the two complex iV— dimensional vectors / and g ac- 
cordingly, 

N 1 r ^ N / N \ 

(/>#) m = S fn9n = (f\g) ► (f>9) () = 2(2 /m/W 9n = {{f\g) 



W n=l \m=l 



In the other words, one may simply abandon the most common definition of the bra 
vector (/| (using the mere operations of the transposition and complex conjugation) 
and replace it by its more general alternative with the symbol ((/| dependent on the 
pre-selected mètric fi — /J > 0. 

The same generalization will also work for Hilbert spaces with N = oo after a 
few necessary but entirely straightforward additional technical qualifications [5, 20]. 
One only has to keep in mind that any modification of the mètric / — > fi specifies a 
different, fi— dependent Hermitian conjugation mapping A — > A*, 

A^^A^, n = vt. (1) 

This means that the class of the corresponding " Hermitian" matrices A = A* varies 
with fi. In order to avoid confusion, let us call the matrices A = A$ quasi- Hermitian 
(= the name introduced in ref. [5]) whenever the positive definite mètrics fi becomes 
different from the trivial identity operator /. 

Any such a specific mètric operator fi > will be denoted by the symbol r\ in 
what follows. We must be careful with this convention since Dirac [22] used the 
same symbol r\ without requiring its positivity. This means that he spoke about a 
pseudo-metric fi, very well exemplified by the indefinite Pauli matrix occurring in 
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the paper by Feshbach and Villars [23] where \x = er 3 ® I. For this reason, authors 
of ref. [5] introduced another symbol for the mètric (77 — > T). Their choice proved 
unfortunate as it causes a highly undesirable confusion with the time reversal symbol 
T so that recently, Ali Mostafazadeh [17, 24] recommended the return to 77 with a 
specific subscript, i] — > r] + . A less clumsy notation was preferred in refs. [25] or 
[14] where one starts from certain auxiliary parity-type indefmite operator V and 
constructs then the necessary alternative, "physical" mètrics as a positive-defmite 
product i] = QV (where Q denotes the so called quasiparity [26]) or r\ = CV (with 
the operator of charge C [18]), respectively. 

2.2 Non-quasi-Hermitian time-evolution generators 

Beyond the "mathematically safe" domain of quantum mechanics as covered by the 
majority of the current textbooks one may find a fairly broad "grey zone" of mod- 
els which do not seem to comply with all the postulates of the theory. One of the 
illustrations may be found in relativistic quantum mechanics [27] where several con- 
sequences of the time-evolution laws seem to be in conflict with the Standard (e.g., 
probabilistic) postulates. For a brief explicit illustration of the possible dimculties 
encountered in similar cases, let us just recollect the exactly solvable Klein-Gordon 
model of ref. [28], 

(td t ) 2 ¥ KG \x,t) =H^¥ KG \x,t), = -d 2 x + m\x) 

where 

9/ . P B 2 — A 2 — Auj B (2A + uj) sinh ux , . 

m 2 (x)=m 2 + - 2 + ^ . (2) 

cosh ujx cosh ujx 

Due to its exact solvability the model possesses the bound-state energies in closed 
form, 



= ±\Jml - (A - nu) 2 , n = 0, 1, . . . , n max , n max = entier[A/u] . 

In the strong-coupling regime with < m < A, a few low-lying states suffer a 
collapse and acquire, formally, complex energies, i.e., Im(E^) ^ 0, im^^) 7^ 
etc. We have certainly left the domain of quantum mechanics. 

Once we stay within the interval of < A < m , the strength of the force remains 
weak and we have the well-behaved eigenvalues and bound states at all the admissible 
indices n. The time-evolution of the system is generated by the "non-Hermitian 
Hamiltonian" of ref. [28]. We may abbreviate ip[ FV \x, t) = i d t ^^ KG \x,t) and 
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(p2 V \x-,t) = ^ KG \x,t) and arrive at the Standard Schròdinger-like time-evolution 
equation 

where the relativistic Feshbach-Villars-type Hamiltonian hS FV ^ is not only manifestly 
non-Hermitian, 

h^f = vh^v-\ v = 

but also non- quasi- Hermitian since detV = — 1. As a "candidate for the mètric" the 
operator V is unacceptable because it is not positively definite. 

In a way coined by Ali Mostafazadeh [29], all the operators A = A* with an in- 
definite fi = V may be called V—pseudo-Hermitian, assuming only that the "pseudo- 
metric" V in our Hilbert-Krein space [30] remains invertible. Keeping in mind the 
paramount importance of the Klein-Gordon equation we probably should accept and 
tolerate the occurrence of the similar time-evolution generators in physics. Recently, 
their use became quite useful and popular in statistical physics [31] and also beyond 
the domain of quantum theory [32]. 




2.3 VT— symmetric Hamiltonians 

Our choice of the symbol V was equally strongly inspired by its role of a pseudo- 
metric and by its coincidence with the specific parity operator in the Bender's and 
Boettcher's pioneering letter [33]. They paid attention to the bound states in complex 
potentials on a line, 

V(x) ~ (ix) 2+s , xeR, Se (-1,2) (3) 

which exhibit an unusual, antilinear parity & time-reversal symmetry, 

VTV(x) =V(x)VT (4) 

where V is parity and the symbol T denotes the time-reflection operator [34]. Under 
the condition (4) the real bound-state energies may be shown to emerge for an 
unexpectedly broad range of parameters [35]. 

Although the symmetry (4) may be found in older literature (cf. refs. [36] and 
[37] working with 5 = 2 and 5 = 1, respectively), the importance of the letter [33] lies 
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in the courage with which its authors opened the discussion of the possible sensibility 
of the study of the non- quasi- Hermitian models in quantum mechanics. They had 
to meet a number of objections a particularly concise sample of which has been later 
formulated by R. F. Streater on his web page [38]. Fortunately, a key breakthrough 
occurred in several subsequent parallel studies of the VT — symmetric Hamiltonians, 
i.e., of the time-generators H such that HVT = VT H or, equivalently, W = VHV 
whenever it is assumed that V^ 1 = V. In these studies people paid attention to the 
exactly solvable differential-equation models with forces (2) [39] or (3) [19] or to 
their linear-algebraic diagonal- matrix reformulations [40]. It has been clarified that 
the "input" pseudo-Hermiticity of a given H may and must be complemented by 
another, positively definite (though, of course, H— dependent) mètric r\. Then, the 
pseudo-Hermitian Hamiltonians re-acquire the quasi-Hermiticity property 

H^rjHrj- 1 , rj = > 0. (5) 

All the requirements are being reduced to an explícit assignment of the second, 
"physical" mètric rj to a given Hamiltonian H. One could speak about an "exòtic" , 
VT — symmetric version of the Standard quantum mechanics. 

3 Point-interaction toy model 

The vast majority of all the applications of VT — symmetric quantum mechanics (PT- 
SQM) will rely upon the assumption that for a given V— pseudo-Hermitian quantum 
Hamiltonian H we find a physical mètric operator 77 which is not too complicated. 
Most often, the consistent theory is expected to be based on the factorized form of 
the mètric 77 = CV defined in terms of an involutive symmetry C = C^ 1 (= "charge" 
[14, 18, 41] or "quasi-parity" [26, 39, 42]). 

The sufficiently simple form of C or i] may only be assigned to a sufficiently 
simple "input" operator H. In this sense, all the VT — symmetric models using point 
interactions [10, 12, 13] represent one of the most natural playgrounds for the explicit 
constructions of r\ or C. In such a context we reported [12] a few purely numerical 
results of the study of a one-dimensional Schròdinger equation 

ip(x) = Eip(x) H±L) = (6) 

with a double-well point-interaction VT — symmetric potential V(x) with complex 
couplings. In our present, non-numerically oriented continuation of this study let 
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us first perforin a trivial re-scaling of the interval of coordinates which replaces the 
usual large though fixed and finite cut-off parameter L ^> 1 by the more comfortable 
value L — 1. Secondly, let us slightly simplify the potential of ref. [12] to a mere 
two-parametric VT— symmetric double-well model 

V(x) = -iÇ6(x + a)+i£6(x-a) xe (-1,1) (7) 

where the potential is reduced to the mere pair of the delta functions proportional 
to a purely imaginary coupling of the size £ and located at a distance measured by 
the variable a G (0, 1). 

3.1 Construction of the exact solutions 

We have to solve eq. (6) with the Dirichlet boundary conditions = and with 

the piecewise-constant potential (7) equivalent to the two matching constraints 

-^ip(±a + 0) - -^ip(±a - 0) = ±i £ i/j(±a) . (8) 

This task will be facilitated by the observation [based on the inspection of eq. (8) 
at small £] that our potential (7) represents a smooth perturbation added to the 
solvable (and safely Hermitian) bound-state problem in an infinite square well which 
possesses the well known and safely real discrete and positive spectrum. This means 
that the proof of the reality and positivity of the energies E n remains trivial in 
the weak-coupling regime. We may assume that £ G (0,£crit), being aware that the 
quantity ^ cr u may vary with both the shift a and excitation n. This makes the explicit 
estimates of the maximal allowed coupling Ç, crü important. Their determination will 
be discussed later. Now, we only put E = k 2 and recollect the usual PTSQM 
normalization convention 

ip(x) ^ i/j s (x)+íi/j a (x) : i/j s (x) = ip* s (x) = ips(-x) , Mx) = ip* A (x) = -ip A (-x) 
which significantly facilitates an ansatz for the wave functions, 

Í>l(x) = (a — i (3) sin k(x + 1), x G (—1, —a) 

ip(x) — < ipc(x) — 7 cos kx + i 5 sin kx, x G (— a, a) (9) 

ÍJrÍx) — (ot + i/3) sin k(—x + 1), x G (a, 1). 
Together with the formula for their derivatives, 

ip' L (x) = n(a — i (3) cos k(x + 1), x G (—1, —a) 

ip'(x) = < ip' c (x) — — K7 smnx + inS coskx, xe(—a,a) (10) 
ip' R (x) = —k (a + i/3) cos k(—x + 1), x G (a, 1) 
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both these formulae enter the matching constraint (8) and transforms it into the 
following four-by-four matrix equation, 

f sin/t(l — a) —cos na \ / a \ 

sin/í(l — a) — sin/ía (3 

— cos/í(1 — a) £ /t" 1 sin k(1 — a) sin/ta 7 

\ £ /t _1 sin k(1 — a) cos/t(l — a) cos na J \5 ) 

It possesses nontrivial solutions if and only if its secular determinant T>{k) vanishes, 



0. (11) 



V(k) = -i jsin 2k + ^ sin 2/ta • sin 2 [/t(l - a)]| = . (12) 

The sequence of the nodal zeros n n = K n (a, £) > of this trigonometric function 
determines the spectrum of the energies E n = > 0. The coefficients a,(3, 7 and 5 
define the wave function and are also given by the nice and compact formulae, say, 

£ sin na sin k{1 — a) . sinrea . . 

7 = ^ " L 6, P = ~. r, s S 13 

k cos k sin /t(l — a) 

f cos na sin k(1 — a) cos na 

8 = -- -7, « = - 7^ r7 (14) 

k sin k sin /t(l — a) 

with an appropriate choice of a convenient normalization. 



3.2 Onset ^ crit of the spontaneous VT— symmetry breakdown 

The compact solvability of our matching conditions enables us to determine the range 
of couplings £ for which all the roots K n = K n (a, £) of our secular eq. (12) remain real. 
We shall see that it happens within fairly large intervals of £ e [0, ^ C rit{o)}i at the up- 
per end of which one encounters the genèric merger and subsequent complexification 
of some of the low-lying energy pairs [30]. 

3.2.1 The reality of the energies in the weak-coupling regime 

In the PTSQM context it is usually emphasized that after the transition through the 
point Çcrit( a ) (which is usually called an "exceptional point" [43]), the complexifica- 
tion of the energies is accompanied by the (so called spontaneous) breakdown of the 
VT— symmetry of the wave functions. The rigorous study of this phenomenon proves 
significantly facilitated in our specific, solvable point-interaction model. Firstly we 
may drop the obvious zero root k — of eq. (12) as spurious, giving just the triv- 
ially vanishing wave functions. Secondly, we may start our considerations from the 



8 



domain of the very small ^ <C 1 where it is easy to observe that the well known 
square-well roots are only slightly perturbed during the growth of £, 



K n (a,0~^ + O(e), n = l,2,.... (15) 

With the further growth of £ we notice that the second term in eq. (12) will be most 
infiuential at the smallest energy levels re 2 . The main uncertainty in an estimate of 
its effect lies in the a— dependence of its trigonometric factor. 

3.2.2 An explícit formula for Ç crit {l/2) 

The mechanism of the spontaneous breakdown of the VT— symmetry will be best 
visible in the most symmetric case with the half-unit displacement a = 1/2 of our 
5— functional interaction points in eq. (7). In the secular equation 



sm re 



í 2 

COS K H (1 — cos re) 

4 re 2 



0, a = \ (16) 



the odd-parity square-well roots [= those with even subscripts n in the zero-order 
eq. (15)] will not move with the growth of £ at all, due to the £— independence of 
the first factor in eq. (16). In contrast, the £— dependence of the odd— n series of 
the roots will be nontrivial and controlled by their implicit definition (16) or, after 
a simplification, by the equation 



u 2 £ 

cos re = — -, /i = -- (17) 

k z — 2 



Such a picture is extremely transparent because the right-hand-side curve is a smooth 
and growing function of re > with a single pole at re = //. Its vàlues are larger than 
one at the smallest re < /i where it cannot produce any roots since we must have 
| cos re| < 1. At all the larger re > /j, — £/2 the latter curve is negative. Of course, 
it should not lie below the minus-one lower limit so that we may set an improved 
estimate of the roots as re > \pl fi = £/\/2. This is a constraint imposed, in effect, 
not only upon the energies re 2 in a given potential but also upon the allowed measure 
£ of the manifest non-Hermiticity at a given energy. 

An inspection of eq. (17) reveals that the growth of £ implies that the first 
two dependent energy levels E\ = re 2 (l/2,£) and E% = reg(l/2,£) move towards 
each other until they coincide at a certain critical strength £ cr j t = 2/io- We have 
/«i (1/2, Çcrit) = re 3 (l/2, Çcrit) — vq > 7r at this point. In the other words, the left- and 
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right-hand-side curves in eq. (17) will touch at a point which is determined by the 
two coupled equations 



COS I/q = 2 2 ' - S111 U = — ^2 



We eliminate 



2 o cos2 ^0 



//„ = -2 i/ 



sm uq 

and get the reduced equation for the remaining unknown, 

v = 2 tan - vq cos z/o- 

We may find its root v m 3.874366817 numerically, at fi m 2.529882472, i.e., at the 
critical strength Ç crit sa 5.059764944. 

The touching-point Ç, cr ü will grow with the excitation. We may conclude that the 
spectrum remains real at all the couplings £ < ^ crit ~ 5.06. This means that the 
PT— symmetry becomes spontaneously broken and some of the energies complexify 
only beyond this critical strength of the non-Hermiticity. 



3.3 Unavoided crossings of the energy levels 
3.3.1 Facilitated construction at a = 1/2 

The touching point of the negative branch of the cosine with an increasing right- 
hand-side hyperbolic curve of eq. (17) is displaced to the right of the first excited 
level, vq > /í 2 (1/2, £) = 7r. The latter level was generated by the sine part of secular 
equation (16) and coincides with the position of the minimum of the cosine curve. 
We may conclude that as long as we had the ground-state level /ti(l/2,0) = n/2 
at £ = 0, the continuous and growing function Ki(l/2, £) of £ crosses the line of 
«2(1/2, £) at some particular coupling strength £ 1(gl2 which lies somewhere inside the 
interval (0,Ç crit ). 

The existence of the similar points of an unavoided crossing of the real energy 
levels has been spotted in several other solvable models [26, 44, 45]. In the present 
a — 1/2 illustrative example these unavoided crossings will involve all the pairs of 
the real levels « 4m _ 3 (l/2, £) and « 4m _ 2 (l/2, £). Their crossings are unique and will 
take place at the special strengths £ = £(4 m -3)®(4m-2) defined by eq. (16) which leads 
to the element ary formula 

£(4m-3)®(4m-2) = V^TT (2m - 1) , 771 = 1, 2, ... . (19) 
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At these crossing points, our non-Hermitian Hamiltonian ceases to be diagonalizable 
and remains only block-diagonalizable, containing the so called Jordan canonical 
irreducible triangular submatrices [46]. Whenever we wish to avoid similar subtleties, 
it suffices to assume that £ < min (£(4 m -3)®(4 m -2)) = £i®2 = \J~2-k ps 4.442882938. 

3.3.2 Approximate constructions in the vicinity of a = 1/2 

The remarkable factorization (16) of our general secular equation hints that in a 
small vicinity of a = 1/2 one could still expect simplifications. Introducing a small 
measure of perturbation a let us write a — 1/2 + <j/{2k) and return once more to 
our exact secular eq. (12) which reads, in the new notation, 

£ 2 

sin 2k H sin (k + a) ■ [1 — cos (k — a)] — . 

In the first two orders of a the second term has the form of the product 



sui h ; d cos h — — a 2 sin k] • (l — cos k — a sin k + - a 2 cos k 



2k 2 V 2 J \ 2 

so that we may study the following approximation of the secular equation 



sin 2k + 



2 



sin k (1 — cos k) — a (1 — cos /t) — - <r 2 sin /t 



C (V) . (20) 



2/t 2 

First thing we notice is that at o ^ the factorizability of eq. (16) is lost so that 
one cannot factor the term sin/t out. This means that all the sine-generated roots 
K2(a,£)? K/ í( a íOí ^(o, £), ••• (which were £— independent at a — 0) may become 
£— dependent even near the square-well regime with very small £ ~ 0. 

Within our precision oí\ + O (a 3 ) let us employ the partially factorized approx- 
imate secular equation (20), 



2 



sin -k { 4 cos -/c cos /t 



2 sin /t sin -k — 2 o sin -k — a cos -/c 
2 2 2 



0. 



2 ( 2 2/t 2 

In its light the subset of the roots /t 4 (a,£), /t 8 (a,£), /íi 2 (a,£), . . . still remains more 
or less £— independent at the sufficiently small a ~ a — 1/2. For the complementary 
family with the safely non-vanishing sin |/t we get, in the first nontrivial approxima- 
tion, the implicit definition 

1 £ 2 

cot -k cos k + - — - [sin K — <j] — 

2 4/t 2 

which indicates that all the roots /t 2 (a,£), n 6 (a, £), /t 10 (a,£), ... which annihilated 
cos(/t/2) at a — 1/2 (i.e., at a = 0) will now grow with o\ As one can easily 
demonstrate, they will all move in proportion to the factor £ 2 /4/t 2 at small £. 
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We may conclude that although the study of our toy model in the domain where 
a — 1/2 does not vanish is perceivably more complicated in the technical sense, a 
straightforward perturbative approach to a systematic derivation of the 0(cr k ) cor- 
rections is still feasible. Hence, we may recommend it as a complement or alternative 
to the less sophisticated numerical constructions. 

Skipping the details we shall now return to a = 1/2 and show how our model 
may be interpreted as a mere small perturbation of its Hermitian square-well £ — > 
limit. 



4 The excitation-dependence of the solutions 

Our construction of the solutions need not necessary be only based on the purely nu- 
merical determination of the roots K n from transcendental equations. Alternatively, 
perturbative methods are also easy to apply, especially when K n = 0(n) is large. 

4.1 A power-series analytic representation of the energies 

Whenever the value of £ stays safely below the VT— symmetry breakdown boundary 
ícrit ~ 5, all the roots K n (a,£) of our secular equation (12) with n — 1, 2, . . . converge 
to their compact square-well limits K n (a, 0) = nn/2 not only when we move £ — > 
but also for the growing n — > oo. In order to describe this phenomenon in more 
detail, let us return back to the secular equation (16) where we may skip the trivial 
case (= even-subscripted closed solutions /t 2m (l/2,£) = mn with m — 1,2, . . .) and 
introduce the following new variables for the odd-subscripted roots, 

(4m-3W Um- 1W 

«4m-3 = ^ l·Xín ' /t4m " 1 = 2 m = 1 ) 2 '···· ( 21 ) 

This simplifies the secular sub-equation (17) which forms the series of equations 

1 Tt (4m-3W + 2x , n 

sm x=— r -, L(m,x) = , m = l,2,..., (22) 

L 2 (m, x) — 1 2/i 

for the unknown x = x m and another, slightly different series for y = y m , 

smy = — -— , C/m,y =^ ^, m = l,2,.... 23 

U z (m, y) — 1 2/x 

The roots remain bounded, a: n < 7r and y m < n, but they depend on both the 
coupling /i = £/2 and index m. As long as equations (22) and (23) are not too 
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dissimilar, we may treat them as two special cases of a single equation, 



sm z 



A/i fi) 



-i 



X = { (4m : 3) " . (24) 



Vm \ + 2 



(4m-l)w 

Once À is small, it is trivial to solve it in closed form, say, by iterations, 

z = AV + A V + 2 AV + ? - A 6 /i 6 + 6 A V + (j + ïlpj A 8 /i 6 + K{\, /i) . (25) 
Here, the O (\ 9 fi 8 ) remainder has transparent structure 

n\ h) = y A V + (36 + AV + (SO + ^pj A V + O (A 13 // ) 

and admits the further higher-order improvements if needed. Note that this series 
is not changing signs which means that the function itself is monotonic and growing 
with A > 0. It behaves like a perturbation series in both A (inversely proportional 
to the excitation number m) and /i — Ç/2 (the coupling) simultaneously. 

4.1.1 A check of consistency 

We may compare the rate of the growth of x m and y m with an increase of the coupling 
/i. This requires a removal of the difference in the m— dependence of the respective 
As, via an obvious transition to the same, shared m— dependent factor r, 

The use of the new variable defines the functions x m = x(t) and y m = y{r) in such 
a way that x(t) = y(—r), i.e., both of them are given by the same series in r, 

4 2 ±4- 3 +Ít4+4V 4 ± (4-2^+241 - 5 +···· 

7T 2 7T 2 ^47T 2 7T 4 ; v 27T 2 7T 5 7T 4 j 

As expected, x(r) grows more rapidly with r so that we arrived at an independent, 
non-trigonometric proof that x rn > y m for all m for which these series converge. 
Unfortunately, in contrast to our previous formula (25) which offers a systematic 
improvement of convergence with respect to both the increasing couplings fi = Ç/2 
and the decreasing excitations m oo, the latter two power series in r have an error 
term 0(fi 2 ) so that they remain less useful unless the coupling fi remains very small. 
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4.2 Matrix elements and wave functions 

As long as we work just with the four-dimensional matrix problem (11), it makes 
sense to derive the perturbed wave functions directly from its a — 1/2 version with 
abbreviations S = sin(/t/2) and C = cos(/t/2), 



í s 





-C 
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a 


\ 





s 





-s 






P 




-c 




s 









7 






C 





c 


) 
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) 



An analysis of such a compactified Schròdinger equation must be performed sepa- 
rately in the two or three different domains of k. 

4.2.1 Special cases of the levels when S = or C = 

Once we admit that S = 0, we infer that a = 7 = and (3 = —5 7^ while 
C = ±1. Wave functions ip(x) ~ sin kx vanish precisely at the points of interaction 
and so they "do not see it" and remain independent of £. The whole family of these 
ip(x) coincides, at all the real non-Hermiticites £, with the pure square-well states 
ift^x) , i/)g(x) , . . . with subscripts n = (mod 4). 

In another special case of the matrix elements such that C = 0, wave functions 
may be re-normalized to having antisymmetric real part, 

Í^l{x) = sin/c(x + 1), xe (-1,-1/2) 

^c(^) = — cos/tx - sin/tx, x e (—1/2, 1/2) (27) 
^fí(x) = sin «(x - 1), xe (1/2, 1). 

At all £ the corresponding energies remain also equal to the Hermitian square-well 
vàlues E 2 ,Eq,.... We may notice that there exists also a non-trivial symmetric 
imaginary part of ip(x) which vanishes for |x| > 1/2 and which grows in proportion 
to the coupling £. 

4.2.2 Genèric Schròdinger equation 

At the genèric matrix elements S 7^ 7^ C our four-dimensional Schròdinger equa- 
tion (26) implies the simplification 5 — (5. This may be inserted in our ansatz (9) for 
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= 0. 



wave functions and enables us to drop the second line in eq. (26) as redundant, 

/ S -C \( a \ 

-C Çk^S S 
k Çk^S 2C j 

This simplifies the derivation of the other two closed formulae, 



P 

V 7 j 



(28) 



[3 = 5 = 



2/t 



tan(/t/2) a , 



7 = tan(/t/2) a . 



(29) 



We are left with the single equation requiring that the secular determinant vanishes. 
We may freely choose a normalization a — 1 and conclude that the determination 
of the wave functions is completed. Their explicit perturbative dependence on the 
excitation number m is mediated by formulae (21) and (25). 



5 The coupling-dependence of the solutions 

Our reduced secular equation (17) contains the trigonometric functions S and C. In 
an opposite direction we may abbreviate g — Ç/k and eliminate S and C, transform- 
ing this secular equation into one of the following two equivalent new formulae, 

1 - Q 2 /2 



C 2 



(30) 



2 - g 2 /2 ' ' 2 - g 2 /2 

These relations re-define the matrix elements in our four-dimensional "perturbed" 
Schròdinger equation (26) as functions of an alternative new variable g which remains 
small in the weak-coupling regime.. 



5.1 An illustration with z = x 



rn 



The choice of z — x m covers the energy levels E\ (= ground state), E5, Eg etc. They 
are all slightly greater than their "unperturbed" , square-well £ = predecessors. 

From our previous analysis of the subset of eigenvalues E4 m _ 3 corresponding to 
z = x m we infer that below the level crossing, i.e., for all the not too large couplings 
£ < £(4m-3)®(4m-2) = V^tt (2m — 1) we may take the square roots in eq. (30) 
unambiguously. Without any loss of generality (i.e., up to an irrelevant overall sign), 
we shall evaluate our matrix elements in accordance with the unique recipe 



S = S x (g) 



2 - g*/2 



> C = C x {g) 



1 - g 2 /2 
\ 2 - g*/2 



> 0. 



(31) 
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In both these functions we may consider the quantity g small and employ their known 
Taylor expansions, 



8 128 1024" 32768 262144 L 
Another Taylor series carrying a weak-coupling perturbation expansion character, 

smcm = i + ^ + 1/ + ^ + ^é>" + o ( e «>) 

enters the formulae (29) which define all the components 7 = 7(0) = 5 and /3 = f3(g) 
of the perturbed wave functions. 

5.1.1 The second special case with z = y m 

The fourth quarter of the set of all eigenvalues comprises the energies E im _i which 
lie below their unperturbed partners and which correspond to the choice of z — y m . 
In this case we have to take the square roots in eq. (30) in accordance with this 
modified expectation, 



S = S v (q) = J jr^nz > -C=-C y {g) = 



2 - g 2 /2 " "~ wW \ 2 - g 2 /2 



l -f' 2 >0. (32) 



The necessary parallel modification of all the previous z = x m results and formulae 
is trivial. 



5.2 Wave functions in dependence on both £ and 1/m 

Our wave functions have to be represented in terms of both the complementary 
and independent measures of the non-Hermiticity À = A(m) and \x = £/2. The 
final transition from the intermediate auxiliary quantity g to this original pair of 
expansion parameters profits from the feasibility of the Taylor-expansion technique. 
Up to corrections O (A 9 /i 8 ) we have the following final result for the matrix elements 
in eq. (26), 

V2 S x [g( A)] - 1 + \f A 2 + \^ A 4 - // A 5 + ^ A 6 - \^ A 7 + (J| f - \ /) A 8 , 
V2 C x [g{ A)] « 1 - \f? A 2 - \^ A 4 + ^ A 5 - H ^ + Z ^ A 7 + (- ±g / + i /) A 8 . 
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This converts our Schròdinger equation (26) of the form Q(X, /i) ip = into its power- 
series four-by-four matrix version 



Q(0, 0) + fi 2 \ 2 Q w + //A 4 Q [2] + ^ 4 A 5 Q [3] + . . .j fi)=0. 

The Standard Rayleigh-Schròdinger perturbation theory [1] might apply and define 
the sequence of corrections, order-by-order in both our small parameters. The exact 
solvability of our problem facilitates the underlying manipulations and reduces the 
abstract algorithm to the explícit formulae (29) where the only ingredient we need 
to know is the expansion of the single auxiliary factor S/C. It is easy to derive its 
sixth-order weak-coupling form 

tan [« m Q, 2^)] » 1 + fi 2 \ 2 + 1 //A 4 - 2 //A 5 + ^ /i 6 A 6 - 8 ^X 7 + (fV ~ ^) A 8 
as well as all its further higher-order improvements whenever needed. 



6 Summary and outlook 

We described an application of the PTSQM formalism to a toy-model Hamiltonian 
H with the real energies E n = k\ which remain non-degenerate for £ < £ 1(g ,2- The 
set of the bound states \n) (9) is not orthogonal since W = V H V ^ H. Still, 
in our Hilbert space we may construct a biorthogonal basis using the second set of 
the eigenstates of corresponding to the same spectrum of course. They will be 
distinguished by the double-bra/ket symbols (e.g., ({n\ or |n)), respectively), being 
most naturally understood as the "left eigenstates of H" [25] defined by equation 

H ] \m)) = E m \m)) = VHV \m)) . 

Due to the involutivity of the parity (V 2 = I) and due to the non-degeneracy of the 
spectrum we know that V\m)) must be proportional to \m) so that we are free to 
define 

\ m )) = Qm^\ m ) 

using any normalization constant q* m (for our future convenience we write here the 
star * which denotes the complex conjugation). It is also easy to deduce the law 
of biorthogonality ((n\m) ~ 5 mn (hint: subtract the two suitably pre-multiplied 
equations H\m) — E m \m) and ((n\ H = E n ((n\). Once we evaluate all the overlaps 
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R n = {(n\n), we may formally write down the explícit decomposition of the unit 
operator, 

OO OO 

1 = E I") TT (H = E l n » H 

n=l - rl ™ n=l ""-n 

as well as the spectral decompositions of our two forms of the Hamiltonian, 

H=E\n)^((nl ^ = J2\n))-^(n\. (33) 

n=l n n n= \ -<V 

In practice, we determine the vàlues of R n by calculating matrix elements of parity 
which are all real but which cannotbe all of the same sign, (n\V\n) = R n /q n . 

The key step now lies in making the PTSQM formalism compatible with the 
postulates of quantum theory by an introduction of some new operator n of "the 
physical" mètric in Hilbert space. We remind the reader that it need not be involutive 
{rf 7^ / in general [5]) but it must be Hermitian and positive definite, i.e., n = 
(in this sense, the indeterminate V is not a mètric). The main purpose of the search 
for such a mètric i] is that it makes our Hamiltonian quasi-Hermitian and, hence, 
physical, 

H ] = r ] Hr ] - 1 . 

Let us now re-interpret this relation as an equation for r\ which we have to solve. In 
the light of eq. (33) we have 

00 E 00 E 

E^N^MH = E \ m ))-^( m \v (34) 

n=l m= i -tt-m 

which implies that, in general, 

00 

V = E \n))M nm ((m\. 

m,n=l 

The backward insertion of such an ansatz in eq. (34) gives, at all n and m, the 
condition E n M nm = M nm E m with the solution M nm = c n 5 nm . The Hermiticity 
and positivity constraints restrict finally the freedom of the choice of the optional 
sequence of c n to the real and positive parameters c n = i] n > 0. Vice versa, any 
choice v^ oyS) of the latter sequence defines an eligible operator of the mètric, 

00 

1N = EI«))^'(H='ÏM>0· (35) 

n=l 

We arrive at a climax and summary of our present message. At any coupling in 
the allowed range £ < £1^2 ~ 4.44 the eigenstates \n) = \n, £) and \n)) = \n,Ç)) 
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pertaining to our toy model H = H(Ç) have been shown to exhibit an asymptotic 
suppression of all the influence of their non-Hermiticity at all the sufhciently large 
excitations, 

|n,£»~M», |n,0~kO>, n>N = N(0- (36) 

The value of the cut-off iV(£) has been shown to grow quite slowly both with the 
non-Hermiticity £ and with the required precision of the approximation (36). We 
may conclude that on an optional level of precision [and assuming that we choose 
Vn"^ = 1 in eq. (35) for all n > iV(£)] we may always treat any pertaining mètric 
V(toy) as an operator which differs from the unit operator just by a finite-dimensional, 
N— parametric separable modification, 

V(toy)*i-Yl IM»«MI+E h£»^ 02/) «^l- (37) 

n=l n=l 

This is the property which our toy model shares with the other non-Hermitian square- 
well models as available and studied in the recent literature [15, 16, 17, 30]. Still, as a 
rule, the majority of properties of the latter models looks much less accessible to any 
non-numerical treatment. In contrast to that, the present energies as well as wave 
functions preserve a compact and closed form, particularly transparent in the most 
symmetric case with a — 1/2. For illustrative purposes this enabled us to construct, 
entirely systematically and in a fully non-numerical manner, an explicit sample of 
the order-by-order corrections in the approximations (36). Thus, our model as well 
as formulae may be expected to find an immediate application in the key separable 
formula (37) which acquires, in this manner, the form of an explicit perturbation- 
series recipe which may be systematically improved up to (in principle, arbitrary) 
precision. 
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